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Abstract 

We re-egzamine the possible magnitude of the super symmetric con- 
tribution to Rf, with imposed all available phenomenological con- 
straints and demanding good quality of the global fit to the precision 
electroweak data. For low tan (3 we find a new region of the param- 
eter space, with M2 ~ |^| and [i < where Rb remains large, 
~ 0.2180 even for the lighter chargino as heavy as 90 — 100 GeV. 
It is an interesting mixture of the up-higgsino and gaugino. The role 
of various phenomenological constraints is discussed in analitic form 
and importance of small but non-negligible left-right mixing in the 
stop sector is emphasized in this context. The large tan/3 option for 
enhancement of Rb is also reviewed. The available data do not rule 
out this scenario. 
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1. INTRODUCTION. 



A considerable excitement has recently been inspired by the Rb and 
R c anomaly 0, ^ ^J, |5|, ]6|, [7|, [|. The succesful tests of the Standard 
Model (SM), to a per mille level ||10|| , are challenged by the measurements 
of the partial widths of the Z° decays into bb and cc quarks which 
disagree with the SM predictions for m t = 180(170) GeV at the level of 
3.7(3.5) and 2.5(2.5) standard deviations, respectively JTTJ. If both results are 
confirmed, the SM and its simplest supersymmetric extension, the Minimal 
Super symmetric Standard Model (MSSM), are ruled out. However, since the 
Rb anomaly is statistically more significant, it is also of interest to discuss 
the possibility of explaining only the larger than in the SM value of Rb. 
Even if R c is fixed to its SM prediction, R c = 0.172 the then measured 
value of Rb = 0.2206 ± 0.0016 is still 3 standard deviations away from its 
SM value. The issue has been addressed in particular in the framework of 
the MSSM [0, |, |, |, |, ^ |, ||. It is well known already for some time that 
in the MSSM there are new contributions to the Z°bb vertex which can 
significantly enhance the value of Rb (but not R c ) if some superpartners 
are sufficiently light (12], |2], |], [5], || ||. More specifically, for low (large) 
tan/3 the dominant contributions are chargino-stop (CP— odd Higgs boson 
and chargino-stop) loops. 

Any improvement in Rb must not destroy the perfect agreement of 
the SM with the other precision LEP measurements and must be consistent 
with several other experimental constraints (which will be listed later on). 
It is, therefore, important to discuss the changes in Rb in the context 
of global fits to the electroweak data (and with all additional constraints 
included). Such fits in the effective low energy MSSM (unconstrained by any 
GUT assumptions about the pattern of soft supersymmetry breaking scalar 
masses) have shown that it is realistic to obtain the values of Rb up to 
R b = 0.2180(0.2190) for small (large) tan (3 values §. Although still away 
by 1.5cr (l.Ocr) even from the central experimental value obtained with R c 
fixed to its SM value, those results provide an interesting improvement over 
the SM prediction R b = 0.2158 (0.2160) for m t = 180 (170) GeV. At the 
same time the overall best \ 2 is smaller than in the SM fits by A% 2 « 4 (5 
for the fit with R c fixed to the SM value) (for m t = 170 GeV). Also, the 
fitted value of a s (Mz) is modified by Aa s (Mz) ~ — 4<5i?b, i.e. lower than 
in the SM fits which give a s (M z ) = 0.123 ± 0.005 0, §, 0, [RJ. This may 
look desirable [0] in view of the apparent hint for some discrepancy between 
the value of a s (Mz) obtained from the SM fits to the precision electroweak 
data and the values a s (M z ) = 0.112 ± 0.005 and a s (M z ) = 0.112 ± 0.007 



obtained from the deep inelastic scattering |L5[ and lattice calculations [|l| 



respectively []. Furthermore, increase of Rb in the MSSM implies some light 
1 The overall average gives a s (Mz) = 0.117 ± 0.005 and includes also the results 
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superpartners, with masses of the order or even below the electroweak scale 
Mz j|, 0, |6], |18| ! (It is may be worth noting small differences between 



various viewpoints: One is to explain the measured value of Rb- The other 
is to take it as a statistical hint for values of Rb somewhat higher then 
predicted by the SM. It is reasonable to take the latter and to consider an 
increase in Rb in the range 0.2170-0.2180(90) as interesting.) 

In this paper, encouraged by the importance of the Rb anomaly for 
experimental search for supersymmetry, we take up this issue once again. We 
clarify certain points of the earlier analysis and provide further insight into 
properties of light superpartners predicted by such global fits. Furthemore, 
we clarify the role of various experimental constraints, including the recent 
new limits on the chargino mass from LEP1.5. Our main new result is that 
for low tan/? the \ 2 of the global fit and the value of Rb depend very 
weakly on the chargino mass (for fixed stop mass) in the range 50—100 GeV. 
The Rb remains at the level of 0.2178 for mc, up to 90 GeV in the 
region where M 2 ~ \fi\ and \i < 0. We also discuss large tan/3 case. 

2. LOW tan/5 REGION. 

Our discussion will be divided into small and large tan (3 cases. To start 
with, let us, however, recall the basic facts from the global fits in the MSSM. 
In order to maintain good agreement of the SM with the bulk of the precision 
data, such as Ap, Mw, sin 2 #f e p{, ... we must avoid new sources of the 
custodial SUv{2) symmetry breaking in the left currents. In the MSSM, 
this is assured when the left squarks of the third generation (and all left 
sleptons) are sufficiently heavy [19|, |J, say, > (9(500 GeV) 0. At the same 
time, an increase in Rb requires a light right stop. So, one needs a hierarchy 

& 

M h » M~ tR or M h » M- t2 (1) 

(in our notation £ 2 denotes the lighter stop) with small left-right mixing. 

The second important fact to recall is the pattern of the chargino sector 
(masses and mixings). The chargino mass matrix 

with 

x= ( „ M 2 V2M w sin(3\ 

\y/2M w cosp n ) 1 ' 

obtained from r decay and jet physics 

2 The actual lower limits on left squarks and sleptons depend on tan/3 value |l9|. f|. 

3 Such a hierarchy is very natural in models where soft scalar mass terms have their 
origin at the GUT scale. Even with universal initial squark mass values to 2 ,, the renor- 
malization group evolution with large top quark Yukawa coupling gives the hierarchy (|l|) 
provided m 2 , » Mf 
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is diagonalized by two unitary matrices Z + and Z T XZ + = 

diag{mc x ,rtic^) with < mc 1 < rric 2 (we follow the convention and 
notation of ref. [[2(J) which determine the projection of the physical two- 
component states \f (i = 1, 2) on the gaugino and higgsino two-component 
weak eigenstates (—iip + , h 2 , —iip~~, ti[) = X~) 

hi = Z^Xt, K = Zl% (4) 



^ = iZl\f (5) 
with the Dirac charginos defined as 



c- = m (6) 




It is important to notice that a physical state may contain different 

admixtures of the "up" and "down" higgsinos (h 2 and hi respectively). 
For instance, it may be almost pure up-higgsino (in which case \Z^\ >> 
\ZV;\) in its lower two-component spinor and almost pure gaugino 

{\Z^\ << \Z^\) in its upper two-component spinor A^- For small values 
of tan/? the pattern of the gaugino-higgsino mixing in the physical states 
is not just determined by the ratio r = M 2 /|//| but also crucially depends 
on the sign of fi. In addition, for fixed M 2 and /i also the pattern of 
the physical masses depends on the sign of /i in a crucial way as shown in 
Fig.l. For better qualitative understanding of those points it is instructive 
to consider the chargino masses and mixing for r = 1. For both signs of the 
parameter and any value of tan (3 the chargino mass matrix @) can be 
easily transformed into symmetric form which is diagonalized by an unitary 
(orthogonal in our case of real [i and M 2 ) matrix i.e. Z' + = Z'_. With 
both eigenvalues positive and ordered, mc x < fnc 2 , (we use the fact that 
tan/3 > 1), the complete diagonalizing matrices read: 
for jj, > 





'-c 6 




Z+ = 


( S0 






\ Sg 


c e ) 




V Cg 


se) 



for jj, < 



( s 9 


-c 6 \ 


Z + = 


( Cg 


-Sg\ 




se J 




\~sg 


-Cg) 



(7) 



(8) 



where s g = sin 9 , eg = cos 9. For the angle 9 we get 
for /i > 

tan 29 = = — < (9) 

^/2M w {sinf3 -cos (3) 
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for n < 



tan 29 = 



2/i 



> 



(10) 



V2M w (sin/3 + cos (3) 



We see that for low values of tan/5 the two cases are very different. For 
/i > we get 9 ~ n/A and cq ~ sq ~ l/y/2. Thus all two-component 
states are mixed. Moreover, the mass eigenstates are 

1 1 M 2 

m Cl2 = fizf M w (sm(3 + cos/3) + -(1 - sin 2/3)— ^ + ... (11) 
v2 4 /i 

for n > My/ (smaller fi values are excluded by experimental constraints) 
and charginos are split in mass by Amc ~ 2Mw- Pure up-higgsino can 
only be obtained for r ^> 1. One can then check that the lighter chargino 
is higgsino-like in both two-component spinors i.e. \Z±\ » \Z±\ and the 
second chargino is heavy. 

For fi < and <^ M w we get 9 ~ and c e ~ 1. The mass 
eigenvalues for low tan (5 are 



First of all, low values of /i are indeed allowed for m Cl > 65 GeV and the 
heavier chargino is still very light. Moreover, the heavier chargino is pure 
up-higgsino (gaugino) in its lower (upper) two-component spinor. For large 
values of r and/or large values of tan/3 the difference between ji > and 
ji < disappears. 

We can study now the supersymmetric contributions to the Z°bb vertex. 
In the low tan /3 region there are two types of relevant diagrams: with stop 
coupled to Z° and with charginos coupled to Z°. Their (renormalized) 
contributions to the vector and axial-vector formfactors of the Z°bb vertex 



F v and F A (C%. = i/> b -f(F v - ^F^Z ) are (in the limit m b = 0) 



given by (the sum is over the two stop and two chargino mass eigenstates): 




(12) 




xfssfiMlM^mc^MiJ 



(13) 




n,m,l 



x/ ss/ (Ml;M fn ,m Cm ,M fi ) 



(14) 
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cp(C) ^ / r}in,l yn.l*yn,m . t- m,Z pn,i* pn,m 

y ~~ 4s w c w ^ -"^cc ^btC^btC + L zcc n btc n uc 

xm Cm m Ct c (m Cm , M^, m C; ) (15) 



^ V - ^ / r -m,Z Tn,l*jn,m , p m >' r>n,l* p 

' 4<Sh/Cw/ 2^ X^zcc^btc^uc + n zcc n b tc n btc 



n,m,l 

r2 



xfffs{Ml,mc m ,M in ,m Cl 



rp(C) ^ ^~" > > / r>m,l jn,l*jn,m j-m,l r>n,l* r>n,m 

~ As w cw i - zcc bic ~ ZCC btC U btC 



As w c w 



xm Cm m Cl c (m Cm ,Mt n ,m Cl ) (16) 

E/ r m,Z yn,l*yn,m p«i,Z -pn,l* pn,r?^^ 
y^zcc^ bicubic ~ n zcc n btc n bic 



n,m,l 

2 



xfffs(M 2 z ,m Cm ,M {n ,m Gl ) 
where svk, cvk stand for sin^vi/, cos 6*^ and the finite functions / ss /, 



and Co are defined in the Appendix B. The couplings read pO 



= r"T ln -^sin 2 ^ n (17) 



L^ c = Z* m *Zj' + 8 ml cos26 w , 

R™£ c = Z^Z 1 ^ + 5 ml cos 29 w (18) 



jn,i _ JL-T ln z 11 - mt T 2n Z 21 

blc s w + V2s w c w M z sin + ' 

R-J = * ^_ T m Z 2U (1Q) 

btc V2s w c w M z cos 



The top squark mixing matrix T*- 7 is defined by: 

t L \_(T xl T 12 \(i l \_(c t -s t \(t x 
tBJ~\T 21 T 22 ){tJ~\s t c t ){t 2 



(20) 



where s t = sin# t , c t = cos 6^. In the case of small tan/3 the couplings 
R™~ c are negligible and we have 5F V = SFa- When the lighter stop, 
i 2 , is dominantly right-handed, as required for a large bi 2 C coupling, its 
coupling to Z° is suppressed (it is proportional to g sin 2 ^)- Therefore, 
diagrams with stops coupled directly to Z° cannot give any significant 
enhancement of Rb- Significant contribution can only come from diagrams 
in which charginos are coupled to Z°. Their actual magnitude depend on 



6 



the interplay of the couplings in the Cft 2 b vertex and the Z^C^Cj vertex. 
The first one is large only for charginos with large up-higgsino component, 
the second - for charginos with large gaugino component in at least one of its 
two-component spinors. As we have seen, this combination never happens 
for fi > 0. Large can then only be achieved at the expense of extremly 
light Cj and t 2 , either already ruled out by the existing mass limits or 
in conflict with global x 2 and/or other constraints such as BR(b — > 57), 
BR(t — > new) . In addition, for fixed mp, and Mf 2 , Rf, is larger for r > 1 
i.e. for higgsino-like chargino as the enhancement of the Cit 2 b coupling is 
more important than of the Z°CiCi coupling. 

For \i < the situation is much more favourable. In the range r « 
1 ± 0.5 the second chargino which for low tan/5, is very close in mass 
to the lighter one (Fig.l), has large up-higgsino component and gaugino 
component. Large couplings in both types of vertices of the diagram with 
charginos coupled to Z° give significant increase in Rf, even for the lighter 
chargino as heavy as 80 — 90 GeV (similar increase in R^ for /i > 
requires mc 1 ~ 50 GeV and Mf 2 « 50 GeV). 

We now turn our attention to a global fit to the precision data and to the 
role of the following constraints: 1) T(Z° — > X1X1) < 4 MeV (in addition 
to the inclusion of the decay mode Z° — > XiXj i n t° the total Z° width 
in the x 2 fit), 2) BR(Z° -> X \xl) < 10~ 4 , 3) M h > 60 (50) GeV for 
small (large) tan/?, 4) BR(b — > 57) in the range (1.2 — 3.4) x 10~ 4 
5) BR(t — > new) < 45% (following ref. pip 0. 6) Recent exclusion 
curves in the (m^, MgJ plane from DO obtained under the assumption 
m Cl > Mi 2 J23|. Since the uncertainty in most of these constraints is 
much larger than in the precision data, we perform the x 2 fit to the latter 
and impose the former as rigid constraints. Importance of the constraints is 
illustrated in Figs. 2 and 3 and should be discussed in more detail. 

The limits T(Z° X1X1) < 4 MeV and BR(Z° -»■ X \xl) < 10~ 4 put a 
constraint on the (M 2 , fi) parameter space provided we make an additional 
assumption about the values of M\ (bino mass). In this paper, for a sake of 
definiteness, we adopt the GUT relation M\ = (5/3) tan 2 8w Then for 
low tan/3 the two constraints eliminate (approximately) a band in the M 2 , 

plane bounded by —50 < < 100 GeV. This is basically due to the fact 
that for M2 and fi in this region neutralinos are too light and/or have too 
strong coupling to Z°. The values of M 2 and fi chosen for Figs. 2,3 
(and 7) are outside the forbidden region. One should remember, however, 
that the limits on Z° — > N®N® and Z° — > N^N 2 are less constraining 
for larger values of M 1; Mi > 0.5M 2 (i.e. heavier LSP). Similarly, the 

4 This is a constraint mainly on the decay t — > t2Nf . In the large tan/3 case, an 
important decay can also be t — > bH + , with M H ± close to Mw- It is difficult 



to distinguish this decay mode from the standard one, t — » bW + |22 and to put an 
experimental upper bound on this branching ratio. 
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bounds on M^ 2 from the DO exclusion curves gradually disappear for 
Mi > O.5M2. Finally, the larger the Mi the better the degeneracy between 
the lighter chargino and the lightest neutralino masses. This is important 
since the new LEP1.5 limit mc, > 65 GeV has been obtained under the 



assuption m Cl — m N o > 10 GeV |24|]. We see that the significance of 



various constraints crucially depends on the ratio Mi/M 2 . With the GUT 
assumption, our results remain on the conservative side. 

The role of the lower limit on the Higgs boson mass (for a compact formula 
for radiatively corrected lighter Higgs boson mass in the limit Ma >> Mz 
see [25]) depends on the mass of the heavier stop and the left-right mixing 
angle. For M^ > 500 GeV (as required for good quality of the global fit) and 
small mixing angles (necessary for large Rb) Mh is above the experimental 
limit Q in a large range of the parameter space. Very small and large left-right 
mixing angles are, however, ruled out by this constraint. This is clearly seen 
in Figs. 2 and 3 where we show the allowed region in the {M$ 2 ,9t) plane 
for fixed M 2 , /i and tan/3. 

The b —>■ sj decay is a very important constraint on the parameters 
space. In addition to the experimental error in the BR(b —>■ 57) (which we 
take at the 2a level) there is large uncertainty in the theoretical prediction 



mainly due to its renormalization scale dependence pq| . In Figs. 2 and 3 
we show the significance of the b — > 57 constraint in two cases: when 
the renormalization scale Q is fixed, Q = nib = 4.7 GeV and with the 
theoretical uncertainty included. In the latter case we consider the result for 
the BR{b — > 57) as acceptable if any of the theoretical values obtained 
with Q varying from nib/ 2 to 2m^ falls into the 2a experimental 
range. Moreover, the results for BR(b — > 57) show weak but nonnegligible 
dependence on the value of a s (Mz). Since the fitted value a s {Mz) depends 
on the change in Rb, in Figs. 2 and 3 for self-consistency we use the value 
a s (Mz) = 0.123 + 5a s where 5a s = —45Rb and the value 0.123 is 
obtained from the SM fit Jl3|, [14|, [Kj. For comparison, in Fig. 2, we 
also show the regions excluded by this constraint with a s (Mz) fixed to 
two different values 0.114 and 0.135. An important message from Figs. 
2 and 3 is the dependence of the b — > S7 constraint on the left-right 
mixing in the stop sector. One should stress that we keep the mass of the 
CP— odd Higgs boson large, Ma = 1 TeV (as needed for large R b ) and 
in consequence the charged Higgs boson is also heavy and its exchange gives 
negligible contribution to the BR{b — > 57). The acceptable values of this 
branching ratio are obtained from the sum \Aw y + Awsy\ of the SM 
W ± exchange and the supersymmetric contribution of the t — C~ loops 
which has to be of the opposite sign. There are two possible solutions: either 



^Important role of the experimental lower bound on Mh in ref. |j] in constraining 
the potential increase of Rb is due to the chosen upper bound < 250 GeV which, 
anyway, looks too low from the point of view of global fit. 
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IAst/sfI ^ l^vT S7 | or I ^susy I ^ l^w^ S7 |- The strong dependence of 
the supersymmetric contribution to the b — > 57 rate on the mixing angle 
$t observed in Fig. 2 and 3 can be understood from the general formulae 
given in [p7f| . For r = l and small angles 9 t (interesting values of ARf, are 
obtained only for small mixing angles 8 t ; large angles are anyway eliminated 
by the constraint from the lighter Higgs boson mass), setting sin 2 9 t ~ 0, 
neglecting the contribution of the heavier stop t\ and using the explicit 
form of the matrices Z + and Z_, eqs. ([7|,|8]), we get the expression: 

A b^ S1 _ n 2„2 M W f (l) 

"^■y/gl, SUSY ~ -Jt 6 2 •> ' 

m Ci 




+ sm29 t s e c e g t 

+ ^f^l^l+'^WAz^] (21) 

where a = sign(fx), g t = m t / \/2Mw sin/3 and the expressions for the 
functions f^] g i{ x ) can be found in [p7|| . Recalling that, for the same values 
of the arguments, functions f^/ g i( x ) are roughly 5 times larger in absolute 
values than the /^/^( x )' s (and both are negative) it is easy to see that for 

/i > Ag^sY is small for 9 t > due to the cancellation of terms of 
order g\ with those of order g t . For small negative 9 t both types of 
terms add up and cancel with the standard contribution making the 

total amplitude unacceptably small. Eventually, for large negative angles, 
the supersymmetric contribution overcomes the standard one making again 
the total amplitude of right magnitude. However this region is excluded 
by the experimental limits on the Higgs boson mass. For negative fi the 
cancelation between terms of order g\ and those of order g t , which makes 
•^susy small enough, occurs for negative angles 9 t . For positive angles the 
supersymmetric contribution is large and cancels the standard one making 
the total amplitude too small. 

Thus, there are two mechanisms for obtaining an acceptable value for 
BR(b — > 57) in the MSSM. One is a cancellation between the H and 



supersymmetric contributions [29]. This mechanism is, however, essentially 
in conflict with the simultaneous increase of Rb due to large negative contri- 
bution of the H ± exchange to Rb- The other mechanism is based on the 
choice of the proper range of the left-right mixing angles. It is certainly inter- 
esting to notice that for small tan (3 the region of acceptable BR(b — > 57) 
partly overlaps with the region of large Rb. Small asymmetry (with respect 
to $t = of the constant AR^ USY contours seen in Figs. 2 and 3 can 
also be traced back to the 9 t dependence of the couplings L^~ c eqn. (TTpj). 
The exclusion curve from the condition BR(t — > t 2 N®) < 45% is also 
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shown in Fig. 2. The dependence of BR{t — > faN®) on the left - right 
mixing in the stop sector, i.e. on the angle 8 t , can be understood from the 
formula for the decay width T(t — > t2Nf): 



s 2 w c^ 64vr 

Jki\2 



1 - 2(x fc + yi) + (x k - y { y 



(K\ 2 + K\ )(1 - **) + 4fte(cf C 



f22l 



where = M~ /m 2 , y^ 
ki i 1 



m 2 N Jm 2 and the couplings read PD] 



Sw^n + c wZp ) T 1 ' c + 



2/ 



life 



4iT-i2fc 



- Sw Z%T 2k + - 
3 sm pMz 



sin (3M Z 

rvAi rplk 

N 1 



Z%T 



(23) 



and Z 3 ^ diagonalizes the neutralino mass matrix p0| . For example, for 
< and r = 1, T(t — ► new) is dominated by T(t — * t 2 N 2 ) because the 
second neutralino has larger up-higgsino component. In that case 



,22 



r 22 _ 



sin (3M Z 
"sin/3M z 



Z^ 2 cos 9 t 



Z N sin ft 



(24) 



and 



m. 



64tt s 2 w c 2 w M 2 sin 2 /3 
x [(1 + 2/2 -^2) - 2 sin 2ft^] 



Z%\ 2 Jl-2(x 2 + y 2 ) + (x 2 -y 2 ) 2 



(25) 



explaining larger BR(t — > new) for negative ft seen in Fig. 2. 

The -DO exclusion contours |23[] in the (M^ 2 , m w o) plane eliminate in 
Fig. 2 a band 70 GeV < M^ 2 < (not shown) and a similar band in 

Fig. 3. This follows from Fig. 10 in ref. pj (in Fig. 2 M 2 = -fi = 58 
GeV, md = 85 GeV and m Nl « 30 GeV (with the GUT assumption)). 
Finally, it is worth recalling that 



*2 



Mf 



I ft 



My 



< 1 



(26) 



where 



M~ = 774 + m 2 + - cos 2/?(M§ - M, 



(27) 



(m 2 ^ is the soft supersymmetry breaking mass term) . If we required m\ > 
then with M h « 1 TeV we get |ft| > 9° for M h ~ 0(100 GeV) and with 
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Mr rs 1.5 TeV — \6 t \ > 6°. The corresponding values of A t are ~ 1 
TeV and ~ 600 GeV respectively. The theoretical status of the condition 
m 2 - > remains, however, unclear and we do not impose it on our global fit 
(with this condition imposed the fit would not have changed too significantly 
in particular given the freedom in the choice of M^ L ). 

The results of a global fit are presented in Figs. 4-6 as projections of x 2 
as a function of for several values of r and two different lower bounds 
for the scan in M^ 2 . We also show the best values of Rb obtained with the 
restriction \ 2 < Xmin + 1 where Xmin * s the minimum of the x 2 f° r fixed 
m Cl . The fitting procedure is the same as in ref. ||. In particular the fitted 
parameters include m t , tan f3, M 2 , fi, a s (M z ), M t - 2 and A t . The values 
of the parameters Mf L , M^o, Mf provided large enough, are not relevant 
for the quality of the fit and we have fixed them at 1 TeV. Similarly, the 
parameters like Mj , My* , masses and mixings of the two first generations 
of squarks to which the fit is not sensitive have been fixed at 1 TeV (in 
the large tan/5 case the best fit is obtained for Mg = 130 GeV due to 
additional contribution of Nf — loops to Rb). 

The value of m t = 170 GeV chosen for the plots is close to the best 
value obtained from the fit. Larger values of m t give worse x 2 an d this is 
a reflection of the well known from the SM fits correlation between the Higgs 
boson and the top quark masses [TU, TJ]. Here the Mh is constrained by 



supersymmetry and cannot follow the increase of m t . In our fit tan/3 is 
bounded so that top quark Yukawa coupling remains perturbative up to the 
GUT scale (we take tan/3 > 1.4 for m t = 170 GeV and tan/? > 1.6 
for m t = 180 GeV). We impose this theoretical constraint to remain on the 
conservative side, as lower values of tan/3 (for the same m t ) give larger 
Rb. For a given m t the best fit is obtained for the lowest value of tan/3 
(for very large tan/3 see next section). One may also wish to impose the 
constraint M 2 > 50 GeV which follows from the lower experimental bound 
on the gluino mass m- g > 150 GeV and the assumed GUT relation for the 
gaugino masses. This constraint does not change the best values of Rb and 
X 2 (but it eliminates e.g. large part of our curve with r = 0.5, fi < 0. 

All discussed earlier constraints are included in those results. The fit 
reveals several interesting facts. For both signs of /x, the value of x 2 
eventually increases towards low values of mc r The reason for this 
behavior of x 2 is the contribution of the neutralinos to Y z which becomes 
too large and correspondingly Uhadr - too small (at this point we should 
again remember about our assumption Mi = (5/3) tan 2 6wM 2 ). This effect 
is much stronger than the decrease of Y z due to the chargino contribution 
to the Z° boson wave function renormalization [|30[| . So, the increase of Rb 
is bounded not only by the lower experimental limits on mc x , M^ but also 
by the quality of the fit to the precision data. Still, for /i > 0, the best fit is 
obtained for mc\ ~ 50 GeV, well below the new experimental limit from 
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LEP1.5, > 65 GeV [2~3|. This limit is valid for mc x — ?njVi > 10 



GeV. In principle, the degeneracy of the chargino and neutralino masses can 
be better than 10 GeV. This occurs for r > 10 with Mi « O.5M2 (note 
however, that supersymmetric contribution to Rb decreases for very large 
values of r) or for any value of r for sufficiently large M\. This is an 
obvious possibility which we do not discuss further. On the other hand, for 
m Cl = 65 GeV the maximal reachable value of Rb is only R b = 0.2173. 

For fi < due to the mechanism explained earlier, \ 2 {Rb) is small 
(large) for much larger values of mc x - In fact the new LEP1.5 limit is 
totally irrelevant in this case. A chargino with mass 70 — 90 GeV and with 
the composition described by M2 ~ |/i| remains an interesting potential 
possibility. The dependence of Rb on the stop mass can be inferred from 
comparison of Fig. 4 and 5. We see that, even for Mj 2 = mp, = 90 GeV, 
Rb > 0.2175. Moreover, a significant enhancement in Rb is consistent 
with both configurations: M^ 2 > mc x and M^ 2 < mc x - In Fig. 1 we 
show the cross sections for chargino production as ^fs = 190 GeV. The 
corresponding cross sections for the production of the pair N^N® are also 
shown. 



3. LARGE tan/?. 

Significant enhacement of Rb is also possible for large tan/3 values, 
tan/3 ~ m t /mb. In this case, in addition to the stop - chargino contribution 
there can be even larger positive contribution from the h°, H° and A 
exchanges in the loops, provided those particles are sufficiently light (in this 
range of tan/3 Mh ~ Ma) and non- negligible sbottom-neutralino loop 
contributions. Recently it has been argued [§] that large tan/3 option 
of an enhancement in Rb is already ruled out by the constraints from 
6 — > ctV t and the absence of 46 events in the Z° decay which should 
have been present as a signature for the bremstrahlung production of a light 
pseudoscalar Higgs boson A at LEP1: Z° — > bbA° — > bbbb. It has been 
claimed that the combination of the two constraints rules out the region in 
the (Ma, tan/3) plane which can give substantial increase of Rb (see Fig. 
2 in ref. ||). It seems, however, that this conclusion cannot be maintained 
after proper account of experimental efficiencies in the search for the process 
Z° bbbb. The decay b -> crV T with BR = 2.69 ± 0.44 % indeed 
puts an upper bound tan/3 < 0.52M#±/(lGe\^) [|3^] which translates into 
a bound on Ma as a function of tan (3 which is given in ref. || . For 
instance, for M4 = 55 (65) GeV, allowed tan/3 < 63 (68). However, with 
proper account of the experimental efficiency in search for 46 events and 
with 10 7 events from LEP1, one cannot expect to get for, say, tan/? = 60, 
a lower limit on Ma from this production mechanism better than about 
40 GeV (since the four 6 tagging procedure is not 100% efficient, the 
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main source of the QCD background are bbgg and bbqq final states) 6 . 
This is well below the reported limit of 55 GeV based on the analysis of 
the process Z° — > A°h° ||33|| . Thus, an enhancement of Rb in the large 
tan p region still remains an open possibility and we briefly summarize the 
relevant results. 

The discussion of large tan (3 case is relatively simpler as most of the 
results is symmetric under the simultaneous change of sign of fi and 9 t . 
The contribution to the Z°bb vertex from charginos-stop loops to the vector 



and axial-vector formfactors are (in the limit m& = 0) given by eqs. ([13 
|T7|) . Compact expressions for the remaining contributions are collected in 
Appendix A. 

In the large tan (3 case, for both signs of \l the chargino composition is 
the same in the up and down Weyl spinors and is a monotonic function of the 
ratio r. For enhancement in Rb the bt^C coupling is more important than 
the Z°CC coupling, so higgsino-like chargino (r ^> 1) is more favourable. 
The chargino masses are then given by 

m C2 - M 2 + + (28) 

The diagonalizing matrices have the form 

— sin 6± cos 9± 
cos 9± sin 9± 

with 



(29) 



2y/2M wf i 2V2M W M 2 
tan2g - = Mf + 2M^-^ tan2 ' + = M|-2M^-^ (30) 

so that indeed for M 2 ^> /i, Mw both angles are close to zero. Notice 
however, that 6 + (whose smallness is important for the coupling fe^C) 
approaches zero much slower than 6L as r is increased. The role of the 
experimental constraints (listed in the previous section) on the parameter 
space is illustrated in Fig. 7 The regions ruled out by the BR{b —> 57) 
are easy to understand. Neglecting the contribution of the heavier chargino 
and stop and restricting to small angles 9 t in the general formula of |27| we 
get for the b — > 57 amplitude the following expression 

Ab->s-y ^ 2 2 a M W Al) ( M h \ 

Ar/* SUSY - ~9t U +riTh/9i \riT) 



6 S. Simion and P. Janot, private communication. Such an experimental analysis has 
not yet been performed. We thank J. Kalinowski, M. Krawczyk, M. Carena and C. Wagner 
for several discussions of this point. 
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• oft ( tanf3 \ ft ft Mw f& ( M l 1 >i- 

l 1 m l A2) ( m 2 t \ 
2< ± J ^ l \M 2 H± ) 

With very light A , the charged Higgs boson contribution is large and adds 
up (with the same sign) to the standard W /± -top exchange amplitude. Thus 
we need a cancellation from the supersymmetric part of eq. ([32|). Due to 
the presence of large tan/? factor, the second term is strongly dominant. 
Therefore, for 9 t < supersymmetric contribution to the b — *■ S7 amplitude 
is of the opposite sign to the standard contribution. Moreover, for angles 6 t 
not very close to zero, the absolute value of the supersymmetric contribution 
exceeds the standard one. Thus, for a given chargino mass there are only 
two very narrow bands seen in Fig. 7 in the plane (8 t , Mi 2 ) where the total 
amplitude is acceptable: in the right one the total amplitude is of the same 
sign as the standard W ± — top and charged Higgs boson amplitudes and in 
the left band it is of the opposite sign. From the Figure 7 it is however clear 
that for large tan f3, the b — > s-y rate does not constrain the value of Rb 
at all §. 

The lower experimental bound on the lightest Higgs boson mass in the 
large tan/5 scenario is ~ 40 GeV. Since in the MSSM M^o « M h o, our 
results are not constrained by this bound. 

In the parameter space which gives enhancement in Rb, also the decays 
t — > new are enhanced. In addition to t — > t 2 Nf , important is also the decay 
t — > bH + . For instance, for m t = 170 GeV, tan/3 = 50, m Cl = 65 GeV 
and M A = 55 (65) GeV we get BR(t -> bH+) ranging from 37(34)% 
up to 49(46)% depending on the stop sector parameters. One should 
remember, however, that the decay t — ► bH + is not easily distinguishable 
from the standard one, t — > bW + , for M H + close to Mw ]22| and we do 
not impose this constraint in the global fit. 

In Figs. 8 and 9 we present the results of a global fit in the large tan (3 
case, together with the corresponding values of Rb. Two important features 
of the global \ 2 is the strong decrease in the quality of the fit for m t = 180 
GeV (compared to the best fit for m t = 170 GeV) and for light charginos 
(below 60 GeV). Those effects are stronger than similar effects for low tan (5. 

The values of Rb are almost insensitive to the value of M^ 2 in the range 
50 — 100 GeV and show the expected dependence on Ma and mc 1 with the 
maximal values for very light charginos. It is, therefore, worth recalling that 
the new limit rrid > 65 GeV is based on the assupmtion mc x — ra^ x > 10 
GeV. As for low tan j3, better degeneracy of the two masses can be achieved 
for r > 10 and/or for Mi > 0.5M 2 . In this case also the quality of the fit 
with charginos below 65 GeV is improved as heavier neutralinos contribute 
less to V z- It is also remarkable, that due to the combined effect of neutral 
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Higgs exchange and the chargino-stop together with the neutralino - sbottom 
contributions, the Rb remains greater than 0.2175 for masses well above 
the present experimental limits. E.g., for m Cl ~ Mf 2 xs Ma ~ 70 GeV, 
R b = 0.2178. 

4. CONCLUSIONS. 

In this paper we re-examined the possible magnitude of the supersym- 
metric contributions to Rb, with imposed all available phenomenological 
constraints and demanding good quality of the global fit to the precision 
electroweak data. For low tan (3 we have found a new region of the param- 
eter space, with M 2 ~ \(J,\ and ji < where Rb remains large, ~ 0.2180, 
even for the lighter chargino as heavy as 90 GeV. It is an interesting mix- 
ture of the up-higgsino and gaugino. The role of various phenomenological 
constraints is disscussed in detail in analitic form and importance of small 
but non-negligible left-right mixing in the stop sector is emphasized in this 
context. 

The large tan (3 option for enhancement of Rb is also summarized, with 
similar conclusions to those presented in the earlier papers. The available 
data do not rule out the possibility of large Rb in the large tan (3 case. 

We conclude that the new LEP1.5 limit, m Cl > 65 GeV still leaves open 
the possibility of a supersymmetric explanation of -Rb up to 0.2180. We 
also conclude that Rb > 0.2175 even for mc 1 = M^ 2 = 90 GeV both 
in small tan/3 and large tan/3 cases. Thus LEP2 may not resolve this 
question. 

Finally, we stress that a good quality of the global fit requires the hier- 
archy M~ t2 < M ix (i.e. M~ tR < M iL ). This hierarchy is natural if the 
low energy values of the soft squark masses have their origin in the renor- 
malization group evolution from the GUT scale with the initial condition 

ml > Mf . 
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APPENDIX A. 



Here we collect formulae for the remaining contributions to the Z°bb 
vertex. 

The contribution from charged Higgs boson is given by: 



5F$ 



v 











[('- 




x b - 












[('- 




x t - 


\ 2 X 



fffs(Ml;m t ,M H ±,m t ) 
m\ co(mt,M H ±,mt) (32) 



5F 



(t) 



8sw c w 



1 - -s 



w 



Xb + -^ S v/Xt 



1 - -s 



w 



Xt + -^s w Xb 



f ffs (M z ;m u M H ±,m t ) 
m 2 t c (m t ,M H ±,m t ) (33) 



(l - 24,) (X b + X,) J,„{Ml; M H ±, m„ M H ±) (34) 



(l - 2s 2 w ) (X b - X t ) f ssf (M 2 z - M H ±,m t , M H ±) (35) 



where we used the abbreviations 



X h 



( m b 



tan (3 



X t 



( m t 



cot (3 



(36) 



\M Z r J ' \M Z 
As is well known, contribution of the charged Higgs boson to R b is negative 



IIJ. Positive contribution to R b is provided by the neutral Higgs bosons. 

3 



SF 



(b,A°) 



V 



tA 3 i l -\ s w) Xb f ffs (M 2 z ;m b ,M A o,m b ) (37) 



5Ff A0) 



X b f ffs (M z ;m b ,M A o,m b ) 



(38) 



SF 



(b,H°,h°) 



i a 3 3 \ o W 



x h 



COS' 



a f ffs (M z ;m b ,M H o,m b ) 



+ sin 2 a /// s (Mf ; m b , M h o,m b ) (39) 
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5F 



(b,H°,h°) 



cos 2 a fff s (Ml;m b ,M H o,m b ) 



+ sin 2 a /// S (M|; m b , M h o,m b ) 



(40) 



(b,A°,H°,h°) 



V 



(41) 



X h 



cos 2 a f ssf (M z ; M A o,m b , M H o) 



+ sin 2 a / SS/ (M|; M A o, m 6 , M,/ 



(42) 



where a is the neutral CP— even Higgs bosons mixing angle |3T 



For completeness we display also formulae for the neutralino - sbottom 
contribution. 



5F^ 



4s w c w ^ 



ml 



T 7-2,771 ( T m 

Zbl V bb 



n*j-l,n , Tjm,n* jjl,n 
bb N bbN bbN bbN 



xf ssf {M 2 z ;M~ b ,m Nn ,Mr 



(43) 



5Ff 



4s w c w ' 



771,2 



t W,m / jm,n* rZ,n t->7ti,71* t->, 



m,n* 7->2,n 



xf ssf (M 2 z ;M im ,m Nn ,M- b[ ] 



(44) 



- 2,m 



->n,Z* T->n,m 



ISwCw 



e 

2s\yCw 



E( r->Z,m t-71,2* t n,m , j t,rn 7->n,t* p7i, 
\ ZNN bb ]\f bb ]\f + ^ZNN n bbN n bbN 
n,m,l 

x fffs(M 2 z ; m Nm , il^ m , m ^ ) 

E/ 7-Z,m 7-n,2* jn,m . pZ,m p7i,Z* pn,m 
\ L ZNN^ bbN ^ bbN + U ZNN n bbN n bbN 

n,m,l 

x m jv m mjy, Co , M s , ) 



(45) 



2swCw 



E(4 



77i j^n,l* j^n,m jj,,m pn,2* pn,m 



n,rn,l 



ZNN bbN bbN 



p7l,t* DTI, 

J ZNN n bbN n bbN 



e 

2s\yCw 



x fffs{M z ; mjv m , M Sm , mjv, 

- 2,777 r 77, Z* r 71,771 



E/ 7-Z,77i 7-n,Z* jn,m p2,777 p7i,2* pn,rra 
yL ZNN L blN lj blN - K ZNN K blN n biN 



n,rn,l 



x mjv m mjv, Co (M|; mjv m , Mr , ) 



(46) 
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where the couplings can be found in PU[ and read: 



L 
R 



l,m 
ZNN 
l,m 
ZNN 



y4l* y4m y3l* y3m 
N N A N N 



y3l y3m* yv y4 

^N^N ^N^N 



741 y4m* 



V, 



l,m 
Zbb 



B ll B lm - -s 2 w S' 



i in 



(47) 
(48) 



- n,l 
J bbN 



R 



bbN 



\/2swCw 



2n r/3l 

N 



B zn z 



V2e 
3cw 



B 2n Z%* + 



^/2s w c w M z cos (3 



In y'il* 
N 



B Z 



(49) 



where the matrix Zn diagonalizes the neutralino mass matrix and the 
matrix B^ for sbottom quarks is defined in the same way as the matrix 
T y ' for stops in eq. (p0|). 
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APPENDIX B. 

Here we give expressions for the functions f ss f and ff ss needed to 
compute SUSY contributions to the Z°bb vertex. 

f ssf (s,m 1 ,m 2 ,m 3 ) = -- + — (a (mi) + a (m 3 ) - 2a (m 2 )) 

s — m? — mo + 2m 2 . , . . . . 

+ -M*.mi^ 3 ) (50) 



2s 

+ \m\ + ^ - m2 ^ m3 - m ^ I c (mi, m 2 , m 3 ) 

~&o(0,m 2 ,mi) - -(mij - mj) b' (m 2 ,m l ) 
^b (0,m 2 ,m 3 ) -^(m 2 2 -ml) b' {m 2 ,m 3 ) 



1 s + mf + m\ - 2m 2 , 
fff s {s;mi,m 2 ,m 3 ) = - H — b [s,mi,m 3 ) 

- — (a (mi) + a (m 3 ) - 2a (m 2 )) 

- W-^W-"!) Co(mi , m2 , m3) (si, 

s 

- -& (0,mi,m 2 ) - -/(m? - m 2 ) 6o(mi,m 2 ) 

- 7&o(0, m 3 ,m 2 ) - i(m 2 - m 2 ) b' (m 3 ,m 2 ) 



Standard two-point functions used in eqs. (|51| , |52|) read: 



2 ' 

m 



a (m)=m ^-l + log— (52) 



Q 



x(x — l)s + xm\ + (1 — x) 



m 



2 



b (s,mi,m 2 ) =r) + I dxlog— ^" ' ' v ^ — Z£lA (53) 

Jo Q z 

Q 2 is the MS renormalization scale and 77 = 2/{d — 4). Derivative at 
s = of &o( s ) mi! %) reads: 

. 1 m? + m 2 , m 2 m\ m\ 

&o(mi ' m2) = + K - mF ° g ^ (54) 

Finally, the c (ml,m2,m3) function is defined as: 

c (ml,m2,m3) = f ^ttt 27777— — 75 2U , 2T (55) 

7 7r 2 [A; 2 — mfJ[(A; + p) 2 — m 2 J[(A; + p + q) 2 — m 3 J 

In the case of the Z°bb vertex (p + q) 2 = Mf and we p 2 = q 2 = m\ « 0. 
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fi (GeV) M (GeV) 

Figure 1: Masses of the two charginos and two lightest neutralinos as a function of /i 
for r — 1 (solid lines) and r = 3 (dotted lines) for tan/3 = 1.4. In the lower panncls 
the corresponding production cross sections of C^Cf and NfN^ at LEP2 (s 1 / 2 = 190 
GeV) are shown. Masses of the sparticles exchanged in the t channel are taken to be 
Mp e = Me a = 500 GeV. For neutralinos also the case Mg H = 50 GeV (dashed and 
dash-dotted lines for r = 1 and r = 3 respectively) is shown. 
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Figure 2: Contours of constant 6R^ USY (solid lines) and various constraints in the 
plane (0 t , MQ for m t = 170 GcV, tan/3 = 1.4, M 2 = -pt = 58 GeV (m Cl = 85 
GeV), M^4 = Mj = 1 TeV. Dashed and dash-doted lines show the b — > S7 constraint 
with different treatement of a s : a s (Rb) denotes the curves obtained with a s (Mz) = 
0.123 — 45 R{j u and with the renormalization scale Q varied in the range (mb/2, 2m&). 
The curve with Q fixed correspond to Q = mt = 4.7 GeV and a s (Mz) = 0.123 (for 
more details see the text). Dotted lines illustrate the Higgs boson mass constraint. The 
allowed region is bounded from below by the BR(t — > new) = 45 % curve and the 
parabolic 6 — > S7 curve a g (i2(,) and from the left- and right- hand sides by the dotted 
curves M/, = 60 GeV. The area below the central dotted curve is also excluded. 
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Figure 3: As in Fig. 2 but for [i > and r = 1.5. [i = 85.5 GeV (m Cl = 55 GeV), 
Ma — Mj = 1 TeV. The allowed region is for t > 0, between the two dashed curves 
denoted by a s (Rb) (with exclusion of the central area 40 < Mh < 60 GeV) and for 
9 t < 0, below the dashed curve a s (Rb) (bounded from the left by the dotted curve 
M h = 60 GeV). 
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Figure 4: \ 2 as a function of mc 1 for r = M 2 /|^| = 0.5 (solid lines), 1 (dashed), 
1.5 (dotted), and 3 (dash-dotted) for both signs of \i for m t = 170 GeV, tan/? =1.4, 
Ma = Mj- = 1 TcV. In lower pannels the best values of with the restriction 

X 2 < Xmin + 1 (here Xmin denotes the best x 2 f° r fixed value of mcj are shown. 
In addition we required Mr > 60 GeV. 



27 



negative ji 



positive id, 




m,= 170 GeV 
tan/3=1.4 
M, >90 GeV 

1 7 



m t =170 GeV 
tan/J=1.4 
M, >90 GeV 



70 80 90 1 00 110 50 60 70 80 90 

m c1 (GeV) m c1 (GeV) 

Figure 5: As in Fig. 4 but with condition Mj 2 > 90 GeV. 
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Figure 6: As in Fig. 4 but for m t = 180. 
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Figure 7: Contours of constant 8R^ USY (solid lines) and allowed by BR(b — ► 57) 
regions in the (8 t , Mj 2 ) plane for m t — 170 GeV, tan/3 = 50, M 2 = 1.5/i, mci = 65 
GeV, M fi = 1 TeV, M A = 55 GeV and M 1r = 130 GeV. Contours of constant 
J2i BR(t — ► t2Nf) are also shown (dotted lines). 
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Figure 8: \ 2 as a function of mc 1 (/U > 0) for r = M 2 /\fJ,\ — 1 (solid lines), 1.5 
(dashed), 3 (dotted), and 5 (dash-dotted) for m t = 170 GeV, tan/3 = 50 and 180 
GeV, tan/3 = 55. M~ ti = 1 TeV, = 55 GeV, M 1r = 130 GeV. In lower pannels 
the best values of i?;, with the restriction x 2 < Xmm + 1 (here Xmin denotes the best 
X 2 for fixed value of mc, ) are shown. In addition we required M^ 2 > 60 GeV. 
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Figure 9: As in Fig. 8 but for M A = 65 GeV. 
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